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SECTION  1 


INTRODUCTION 

This  report  completes  the  technical  effort  by  the  General  Electric  Space 
Division  for  the  Air  Force  Geophysics  Laboratory  relative  to  the  Design  Study  of 
a  Laser  Radar  System  for  Spaceflight  Application.  The  primary  objective  of  this 
study  Is  to  Investigate  conceptually  the  requirements  for,  and  formulate  a 
realistic  design  for,  a  laser  radar  system  which  Is  compatible  with  flight 
operation  aboard  the  Space  Shuttle  System  or  other  AF  satellite,  including  also, 
a  proof -of -concept  balloon-based  experiment.  The  primary  scientific  purpose  of 
the  spaceflight  experiment  is  the  determination,  on  a  near  real-time  basis,  of 
the  atmospheric  density  at  all  altitudes  of  interest,  i.e.,  40  km  to  the  ground, 
while  at  the  same  time  retaining  as  a  goal  a  minimum  of  hardware  complexity. 

The  primary  objective  of  this  study  has  been  accomplished,  in  that  an  operational 
concept  of  a  laser  radar  system  for  spaceflight  applications  has  been  developed. 

The  baseline  concept  uses  the  fundamental  and  frequency-tripled  neodymium 
wavelengths  in  combination  with  the  difference  in  wavelength  dependence  of  Rayleigh 
scattering  (i/A4)  and  Mie  scattering  (  +*»  1/X  )  to  separate  the  neutral  atmosphere 
returns  from  those  of  aerosols  and  particulates.  A  summary  of  this  design  study 
which  led  to  the  selection  of  the  baseline  lidar  system  design  concept  is  provided 
in  Section  2  of  this  report. 
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Section  3  provides  a  description  of  the  study  objectives  with  Section  4  pro¬ 
viding  the  approach  to  ensure  that  these  objectives  can  be  met  with  a  viable 
instrument  design  concept. 

Shuttle/Spacelab  experiment  performance  results  are  provided  in  Section  5  for 
both  the  Rayleigh/Mie  and  Differential  Absorption  Lidar  (DIAL)  analysis  techniques 
studied  for  the  determination  of  atmospheric  density. 

Section  6  provides  the  balloon  experiment  performance  results  based  on  the 
Rayleigh/Mie  analysis. 

Experiment  hardware  descriptions  are  provided  in  Section  7  for  the  implementation 
of  an  operational  atmospheric  density  measurement  concept  within  the  Atmospheric 
Lidar  Multi-User  Instrument  System  (ALMIS),  the  Standard  Test  Rack  (STR)  Lidar 
Experiment,  as  an  adjunct  to  the  WINDSAT  Lidar  Experiment,  and  as  a  proof-of- 
concept  (feasibility  demonstration)  Balloon-Based  Lidar  Experiment. 

This  report  is  published  in  a  facing  page  format  in  which  a  summarizing  or 
pictorial  presentation  is  provided  on  the  righthand  page,  while  the  facing 
lefthand  page  contains  a  detailed  narrative  of  the  points  summarized  and/or 
illustrated  on  the  facing  page.  This  format  is  found  most  useful  to  facilitate 
the  use  of  this  material  in  technical  and  management  presentations. 
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APPROACH 


Our  approach  to  this  study  is  two-fold:  first,  to  consider  and  evaluate 
several  density  measurement  techniques  which  will  provide  the  desired 
accuracy  and,  second,  to  develop  potential  experiment  configurations 
which  can  be  used  for  implementation  of  the  desired  technique(s) .  As 
a  result,  the  following  portions  of  this  section  are  divided  in  two 
parts  with  the  first  portion  providing  our  approach  towards  satisfying 
scientific  requirements,  and  the  second  portion  providing  our  approach 
towards  satisfying  systems  engineering  requirements,  stressing  in 
either  case,  an  approach  which  will  provide  a  viable  operational 
experiment  concept  for  spaceflight  application. 
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PREFERRED  CONCEPT 
(FEASIBLE  AND  COST  EFFECTIVE) 
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SHUTILE/SPACKLAB  EPBRDffiWT  PgRFOgMAgCE 


This  section  is  divided  into  the  following  three  major  parts: 


1)  Computer  Simulation  Program  -  this  subsection  presents  the  key  elements 
used  in  the  computer  simulations  for  the  analysis  of  the  Raylelgh/Mie 
and  02“DIAL  density  measurements  techniques. 

2)  Rayleigh/Mie  Analysis  -  this  subsection  presents  the  equations  and  results 
obtained  to  evaluate  the  expected  performance  of  this  technique  from  the 
Space  Shuttle. 

3)  02-DIAL  Analysis  -  this  subsection  presents  equations  and  expected 
accuracy  results  as  they  pertain  to  this  measurement  technique. 
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ATTENUATION  COEFFICIENTS  FOR  AEROSOL  TRANSMITTANCE 
(ABSORPTION  AND  TOTAL  EXTINCTION) 


®  ATTENUATION  COEFFICIENTS  FOR  AEROSOL 

TRANSMITTANCE 

GENERAL  (ABSORPTION  AND  TOTAL  EXTINCTION) 

ELECTRIC 


2.0  3.0  4.0 

WAVELENGTH  (pm) 


SIZE  DISTRIBUTIONS  FOR  LOWTRAN  3B  AEROSOL  MODELS 


SIZE  DISTRIBUTIONS  FOR 
LOWTRAN  3B  AEROSOL  MODELS 


aerosol  models  comparison 
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; 


WAVELENGTH  (microns) 


URBAN  AEROSOL  MODEL 


WAVELENGTH  (microns) 


NORMALIZED  SCATTERING  PHASE  FUNCTION  OF 
AEROSOL  AND  AIR  BASED  ON  THE  MODEL 


(1971),  Optical  Properties  of  the  Atmosphere  (Revised) ,  AFCRL-71-0279. 


)  60  60  100  120  140  160  ISO 

SCATTERING  ANGLE  4  (d«g) 


s MATTERING  angle  PHASE  FUNCTION  FOR.!. 


SCATTERING  ANGLE  PHASE  FUNCTION 
FOR  1.06  MICRONS 


.1 ,jl: 


®"5 

*  Ml 
Ml  ^ 


CD  Ml 
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SCATTERING  ANGLE  9 


SCATTERING  ANGLE  PHASE  FUNCTION  FOR  10.6  m 


This  figure  identifies  the  aerosol  properties  for  operation  at  10.59 
micrometers.  These  data  are  included  for  reference  purposes  only. 
Our  model  is  capable  of  handling  backscattering  from  the  CO2  laser 
at  lO.b  microns,  should  it  be  desired  in  some  future  analysis. 
However,  we  are  not  considering  operating  with  a  CO2  laser  at  this 
t  ime. . 


*  The  reference  source  is  on  page  72. 
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ATMOSPHERIC  AEROSOL  MODELS 


This  figure*  provides  the  vertical  distribution  of  the  attenuation  coeffici¬ 
ents  for  various  aerosol  models  which  have  been  developed  for  the  boundary 
layer,  the  upper  troposphere,  the  stratosphere,  and  the  upper  atmosphere. 

A  good  description  of  these  models,  and  some  examples  of  scaling  to  other  wave¬ 
lengths  can  be  found  in  the  reference  source  for  this  figure  from  which  the 
following  summary  has  been  extracted. 

(.a)  For  the  Boundary  Layer  (below  2  km),  10  models  have  been  defined 
which  describe  the  aerosols  in  rural,  urban,  and  maritime  environments  for 
several  surface  meteorological  ranges  between  2  and  50  km. 

(b)  For  the  upper  troposphere  there  are  two  models  which  represent  spring 
and  summer  conditions  versus  fall  and  winter  conditions. 

(c)  In  the  stratosphere  (up  to  30  km),  models  are  presented  for  background, 
moderate,  high,  and  extreme  volcanic  conditions  for  each  of  the  two  seasonal 
models. 

(d)  For  the  upper  atmosphere  (above  30  km),  two  models  are  presented. 

One  of  these  corresponds  to  the  most  likely  background  conditions,  and  the 
other  represents  the  high  aerosol  concentrations  often  observed  at  these 
altitudes  (in  thin  layers) . 

Also  shown  for  comparison  are  the  Rayleigh  profile,  and  ELTERMAN's  1968  model. 

Atmospheric  Aerosol  Models  used  in  our  analysis  are  as  follows: 

(a)  Temperature/Pressure 

-  Mid-latitude  summer 

(b)  Aerosol  Concentration 

-  Standard  Rural 

-  Visibility 

Variable 
23  KM 

*  The  reference  source  is  on  page  72. 
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ALTITUDE  (Km) 


ATMOSPHERIC  AEROSOL  MODELS 


NUMERICAL  TECHNIQUE?;  USED  IN  OUR  ANALYSIS 


The  numerical  techniques  summarized  below  were  used  for  the  design  of  the 
laser  radar  systems  considered  in  this  study  and  are  based  on  the  scattering 
models  implemented  in  the  LOWTRAN  3B  computer  program.  Atmospheric  trans¬ 
mission  is  computed  by  integration  of  the  atmospheric  extinction  coefficient 
along  the  laser  beam.  Absorption  and  scattering  coefficients  are  derived 
from  the  tabular  atmospheric  properties  provided  by  LOWTRAN.  Interpolation 
is  performed  linearly  in  temperature  and  logarithmically  in  pressure  and 
species  concentration.  The  integration  assumes  that  the  extinction  coefficient 
can  be  fit  by  a  quadratic  in  pressure  and  that  the  pressure  decays  exponentially 
with  altitude  over  an  integration  step.  The  step  sizes  used  for  the 
calculations  are  1  km  at  low  altitudes  and  2  km  for  altitudes  above  35  km. 

The  top  of  the  atmosphere  is  assumed  to  be  86  km. 

For  the  0 DIAL  calculations  the  details  of  the  absorption  line  profile 
are  required.  These  were  obtained  by  assuming  Voight  line  shapes  with  a 
spectral  resolution  of  0.001  cm  1  at  the  line  center.  The  resolution  was 
doubled  at  every  other  wavenumber  point  and  the  line  wings  were  cut  off 
at  10  cm  The  laser  beam  was  convolved  with  the  absorption  line  by  inte¬ 
gration  using  Simpson's  rule. 


•  ALTITUDE  RESOLUTION 

EXPONENTIAL  INTEGRATION 
1  -  2  KM  BITES 
65  POINTS  IN  86  KM 

•  WAVENUMBER  RESOLUTION 

.  LINE  CENTER  0.001  CM-1 

.  SPACING  DOUBLES  EVERY  OTHER  POINT  (SIMPSON'S  RULE) 
LIMIT  2  ABSORPTION  LINES 

•  CONVOLUTION  WITH  LASER  PULSE 

SIMPSON'S  RULE 
LINEAR  AT  EDGE  OF  PULSE 
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RAYLEIGH /MIE  ANALYSIS 


LIDAR  EQUATION 


86 


The  following  pages  outline  the  application  of  the  LIDAR  equation  to  the  measurement 
of  atmospheric  density  by  the  Rayleigh/Mie  separation  technique. 


LIDAR  EQUATION 


EXTINCTION  COEFFICIENT  EXCLUD¬ 
ING  ABSORBING  GAS 


ANALYSIS  OF  RAYLEIGH/MIE  SEPARATION 
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PHOTON  RETURNS  AT  15  KM  USING  LOWTRAN  3B  ATMOSPHERE 


This  figure  Illustrates  the  potential  strength  of  the  Rayleigh/Mie 
approach.  At  15  km  the  return  signal  at  1060  nm  Is  dominated  by  Mie 
(aerosol)  scattering.  Rayleigh  (molecular)  scattering  contributes  only 
25%  of  the  total  return.  Thus,  with  even  a  rough  knowledge  of  molecular 
concentration  (density),  we  can  measure  aerosols  very  accurately. 

At  353  nm,  the  aerosol  effects  are  strongly  overshadowed  by  molecular 
scattering,  and  contribute  only  10%  of  the  total  signal.  This  contribution 
is  corrected  for  by  the  1060  nm  aerosol  measurement  to  give  accurate  density 
values . 

A  value  of  .04  was  used  for  the  normalized  aerosol  backscatter  phase  function 
0A.  This  value  is  used  in  the  following  series  of  figures  for  the  Rayleigh/Mie 
analysis  of  the  Shuttle/Spacelab  experiment. 


EXPECTED  LIDAR  RETURNS  FOR  STANDARD 
LOWTRAN  3B  ATMOSPHERE  (X  "  353  NM) 


'.■Then  all  the  atmosphere  and  aerosol  parameters  are  combined  with  a  Shuttle-borne 
Lidar  system,  the  computer  model  predicts  the  353  nm  signal  levels  shown  In  the 
accompanying  figure. 

The  "Total"  return  is  the  actual  count  that  would  be  recorded  in  telemetry.  The 
contributions  from  "Aerosol"  and  "Molecule"  are  defined  by  the  ratio  of  aerosol 
and  Rayleigh  backscatterig  coefficients  at  each  altitude. 

The  curves  displays  several  interesting  features; 

A)  Total  signals  per  shot  are  not  exceedingly  high.  It  is  clear  that 
multiple  shot  averaging  will  be  required  to  achieve  desired  accuracies. 

B)  Above  20  km,  aerosols  contribute  less  than  10%  of  the  total  signal  at 
353  nm.  Therefore,  even  fairly  rough  aerosol  measurements  at  1063  nm 
would  be  sufficient  to  achieve  10%  density  measurements  above  20  km, 

if  enough  counts  can  be  measured  to  provide  a  statistically  significant 
signal . 

C)  The  stratospheric  aerosol  layer  in  the  7-20  km  region  shows  up  clearly 
in  the  simulation.  The  range  of  variability  of  this  layer  was  shown 
earlier  in  the  figure  on  page  83.  Uncertainties  in  the  properties 

of  these  aerosols  is  a  significant  potential  error  term  in  the  measure¬ 
ment  . 

D)  The  boundary  layer  aerosols  near  the  ground  create  the  most  dramatic 
feature  in  the  simulation,  as  the  aerosol  return  overwhelms  the  molecular 
return  below  5  km.  Also  in  this  region,  total  extinction  effects 
absorb  and  scatter  up  to  85%  of  the  return  signal.  The  combination  of 
increased  aerosol  return  and  significant  total  extinction  severely 

limit  the  achieveable  accuracies  below  5  km. 
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EXPECTED  LIDAR  RETURNS  FOR  STANDARD 
LOWTRAN  3B  ATMOSPHERE  (Aa  1060  MM) 


The  1060  signal  provides  a  good  measurement  of  aerosols  over  the  entire 
altitude  range.  Signal  levels  are  useable,  and  aerosols  contribute  the 
majority  of  the  signal  that  is  observed.  The  influence  of  the  ground 
layer  is  again  significant,  even  in  a  clear  model  with  23  km  visibility 
at  the  ground.  Both  aerosol  scattering  and  total  extinction  effects 
increase  sharply  below  5  km  ,  varying  by  factors  of  two  or  more  over  a 
single  1  km  range  bin.  This  will  make  useable  data  extremely  difficult 
to  extract  reliably  from  the  actual  Lidar  return. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
(SINGLE  SHOT) 


The  expected  density  uncertainties  from  a  single-shot  Raylelgh/Mie 
measurement  from  Shuttle  are  shown  in  the  accompanying  figure. 

Shot  noise  at  353  nm  is  the  dominant  source  of  error  at  high  altitudes. 

At  mid-altitudes,  shot  noise  and  ratio  uncertainties  are  approximately  equal 
contributors.  Below  5  km,  the  increase  in  aerosol  scattering  cause  rapid 
increases  in  several  error  terms,  with 


growing  from 


^  ( 21 )2 


or  lower  up  to  (502)  near  the  ground. 


Similarly 


near  the  ground.  These  two  effects  appear  to  provide  fundamental  experiment 
limitations  at  low  altitudes. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
(100  SHOTS) 


When  100  shots  are  averaged,  the  performance  of  the  Lidar  Increases 

significantly.  With  a  reasonable  understanding  of  aerosol  scattering  ratio, 

(  S  Ratio/  )  a*  50%,  we  can  achieve  better  than  10%  density  accuracy 
/  Ratio 

from  5  km  to  above  30  km. 

With  our  baseline  laser  firing  at  10  pulses  per  second,  100  shots  takes 
10  seconds.  Typical  orbit  velocity  is  7km/ sec,  so  we  must  average  over 
a  70  km  ground  trace. 

Our  linearized  error  analysis  is  not  really  comprehensive  enough  to  address 

residual  errors  in  the  2%  to  5%  range.  However,  we  are  confident  that 

a  more  extensive  iterative  analysis  would  be  able  to  demonstrate  1-2% 

errors  over  most  of  the  upper  altitude  range  for  (  ^Ratio/  small. 

/  Ratio 
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DIFFERENTIAL  ABSORPTION  LIDAR 
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return.  The  error  for  off-line  returns  for  a  given  range  cell  is  then  calculated  by  Poisson 
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EXPECTED  LIDAR  RETURNS  -  02  DIAL 
(ON-LINE  FREQUENCY  CENTERED  ON  SADDLE) 


The  02  DIAL  experiment  was  simulated  on  our  computer  using  the  LONTRAN 
3B  model  atmosphere.  The  saddle  region  between  the  R  Branch  pair  absorp¬ 
tion  lines  (J"  •  13  and  J"  =  15)  was  analyzed  relative  to  the  laser 
line  widths  shown  by  the  curve  data  in  the  accompanying  chart. 

The  "OFF-LINE"  curve  provides  the  nominal  photon  returns  from  broadband 
aerosol  and  molecular  scattering  as  a  function  of  the  range  bin  altitude. 

The  three  (3)  "ON-LINE"  curves  show  the  effects  of  0^  absorption  on  the 
photon  return  signals  for  a  laser  line  frequency  centered  in  the  saddle 
region, also  as  a  function  of  the  range  bin  altitude.  "ON-LINE"  returns 
were  calculated  for  laser  line  widths  of  0.1,  0.5,  and  0.75  cm  ^  (rectangular 
line  profiles  were  used) .  The  reduced  photon  returns  from  broader  lines 
are  commensurate  with  the  width  of  the  saddle  reg  on  (e.g.,  see  figure 
on  page  123  which  shows  the  saddle  width  for  the  f  km  altitude  bin). 

The  large  attenuation  of  the  "ON-LINE"  return  signals  at  low  altitudes 
is  naturally  due  to  a  shrinking  saddle  amplitude  with  decreasing  altitude; 
while  very  little  0^  absorption  is  seen  above  15  km,  due  to  the  combination 
of  a  shrinking  0^  absorption  line  width,  as  altitude  increases,  and  a 
saddle  amplitude  equal  to  the  "OFF-LINE"  amplitude. 

The  "ON-LINE"  return  signals  for  a  laser  line  centered  on  the  saddle  are 
strong  enough,  and  sufficiently  separated  from  the  "OFF-LINE"  returns, 
to  be  promising  at  the  lower  altitudes  (  ^  10  km) ,  even  though  they  are 
somewhat  sensitive  to  variations  in  the  laser  line  width.  However,  at 
the  higher  altitudes  (  > 10  km),  differences  between  "ON-LINE"  and  "OFF¬ 
LINE"  absorption  (AT)  become  indistinguishable,  and  therefore,  meaning¬ 
ful  density  determination  measurements  are  essentially  unobtainable. 
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EXPECTED  LIDAR  RETURNS  -  C>2  DIAL 
(ON-LINE  FREQUENCY  CENTERED  ON  PEAK  ABSORPTION  LINE) 


Simulations  were  also  run  for  returns  from  a  laser  line  centered  on  the 
J"  =  14  line  at  7600.4 93A,  in  order  to  explore  the  possibility  of  working 
on  a  peak  absorption  line  instead  of  in  the  saddle  region. 

The  figure  on  the  accompanying  page  shows  the  "ON-LINE"  photon  returns  versus 
altitude  for  a  range  of  laser  line  widths  of  0.1,  0.15,  0.2,  0.25,  0.5,  and 
0.75  cm  ,  including  the  curve  of  the  "OFF-LINE"  returns,  which  is  the  same 
curve  as  shown  in  the  preceding  chart. 

Since  the  center  of  the  peak  absorption  line  is  totally  black,  even  at 
30  km,  the  "ON-LINE"  curve  data  provides  the  total  absorption  of  a  line  with 
a  finite  width,  rather  than  the  line-center  absorption,  which  would  have  no 
return  at  all.  As  a  result,  increasing  the  laser  line  width  will  increase 
the  return  signal,  since  greater  portions  of  the  saddle  and  OFF-LINE" 
regions  will  be  included  in  the  returns;  as  a  consequence,  the  absorption 
information  in  the  measurement  will  be  diluted.  This  is  particularly  true 
at  high  altitudes  where  the  width  of  the  peak  absorption  line  is  very 
narrow  (the  figure  on  page  125  serves  as  a  typical  example). 

Although  narrower  "ON-LINE"  laser  line  widths  will  work  better  in  the  high 
altitude  regions,  their  utility  rapidly  diminishes  in  the  low  altitude 
regions.  This  is  readily  seen  it'  the  "ON-LINE"  curve  data  where  the  photon 
returns  rom  narrower  line  widths  are  extremely  sensitive  in  the  low 
altitude  regions. 

A  predominant  conclusion  which  can  be  drawn  from  the  data  is  that  the  laser 
must  be  capable  of  optimizing  the  line  width  for  each  altitude  of  interest 
if  meaningful  measurements  are  to  be  obtained  throughout  the  range  of 
altitudes  shown  in  the  chart.  The  impact  of  laser  line  widths  is  brought 
out  in  the  error  analyses  provided  in  the  remainder  of  this  section. 
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EXPECTED  ERRORS  -  02  DIAL 
(ON-LINE  FREQUENCY  CENTERED  ON  SADDLE) 

Using  the  error  equation  shown  on  the  preceding  page  with  the  curve 

data  on  page  127,  errors  were  developed  for  "saddle-region"  02  DIAL 

measurements  from  Shuttle,  and  are  shown  by  the  curve  data  in  the 

accompanying  figure.  It  should  be  noted  that: 

o  Single  shot  means  one  on-line  pulse  and  one  off-line  pulse 

from  a  tunable  dye  laser,  i.e.,  an  output  of  1J,  in  each  pulse. 

o  Laser  line  width  is  modelled  as  a  rectangular  line  with  constant 
total  energy.  The  three  widths  bracket  the  probable  optimum 
width  based  on  the  absorption  profiles  described  previously. 

The  curve  data  portrays  three  general  regions: 

1.  High  altitudes  (  >  12  km)  where  oxygen  absorption  in  the 
saddle  is  so  weak  that  the  AV  term  in  the  error  equation 
dominates. 

2.  Middle  altitudes  (2-12  km)  where  absorption  and  return 
signal  are  optimally  balanced. 

3.  Low  altitudes  (^2  km)  where  strong  absorption  reduces 
return  signal  strength  to  levels  where  1/P  shot  noise  begins 
to  dominate  error. 

In  addition,  the  data  clearly  shows  that  this  concept  (02  saddle- 
region  DIAL)  is  extremely  sensitive  for  obtaining  accurate  density  measure¬ 
ments  in  the  10-30  km  altitude  region,  which  is  a  major  goal  of  the 
Shuttle/Spacelab  LIDAR  experiment,  i.e.,  even  with  100  shot  averaging, 
the  expected  errors  would  be  greater  than  the  desired  measurement  accuracy 
goal  of  10%.  However,  this  technique  looks  very  promising  for  altitude 
regions  from  10  km  to  ground  level  and,  therefore,  can  be  combined  with 
the  Ravleigh/Mie  technique  (see  page  109)  to  provide  measurement  accuracies 
of  beti  er  than  10%  (100  shot  data)  for  all  altitudes  of  interest,  i.e., 
from  the  ground  to  30  km. 


EXPECTED  ERRORS  -  02  DIAL 

(ON-LINE  FREQUENCY  CENTERED  ON  PEAK  ABSORPTION  LINE) 


In  order  to  achieve  stronger  absorption  at  higher  altitudes,  and  re¬ 
duce  the  1/t/r  contribution  to  density  error,  an  assessment  was  made 
of  the  errors  obtained  when  operating  the  "ON-LINE"  frequency  directly 
on  the  center  of  a  peak  absorption  line  (J"  =  14)  of  the  "saddle"  pair. 

As  shown  in  the  figure  on  the  accompanying  page,  a  broad  range  of  laser 
linewidths  was  evaluated  across  the  altitude  region  of  interest. 

The  very  narrow  line  (0.1  cm  ^)  is  so  strongly  absorbed  by  upper 
atmospheric  0  that  returning  signals  decrease  sharply  with  altitude 
(increasing  optical  thickness).  The  broadest  line  (0.75  cm  )  loses 
little  energy  to  absorption  at  high  altitude,  leading  to  large  1/A^ 
errors.  Intermediate  linewidths  (0.15  -  0.25  cm  ^)  display  error  pro¬ 
files  optimized  for  the  10-20  km  region,  but  are  too  strongly  absorbed 
to  provide  low  altitude  returns. 

The  curve  data  in  t  is  figure  (and  also  in  the  preceding  figure)  show  that 
the  02  DIAL  technique  can  provide  density  measurements  within  the  10% 
accuracy  goal  (multiple  shot  averging)  over  specific  altitude  regions. 
However,  a  laser  system  must  be  developed  with  the  capability  to  cover  a 
range  of  tunable  line  frequencies  and  linewidths  in  order  to  provide  data 
over  the  total  altitude  region  of  interest.  Such  a  system,  ruggedized  for 
an  unattended  space  mission,  is  currently  not  available  for  a  near-term 
application. 
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EXPECTED  ERROR  SENSITIVITY  TO  LINEWIDTH  STABILITY 


(ON-LINE  FREQUENCY  CENTERED  ON  PEAK  ABSORPTION  LINE) 


The  sensitivity  of  density  error  to  laser  linewidth  stability  is  shown  in  the 
figure  on  the  accompanying  page.  The  effects  of  minor  variations  in  linewidth  are 
referenced  to  a  nominal  0.2  cm  ^  laser  linewidth  corresponding  to  the  data  shown  in 
the  preceding  figure.  Density  errors  were  analyzed  for  1,  2,  5  and  10%  linewidth 
variations,  assuming  no  1/P  shot  noise  (i.e.  multiple  shots).  These  errors  should 
be  RSS  with  those  on  the  preceding  figure  if  combined  errors  are  to  be  obtained. 

For  the  cases  considered,  variations  in  laser  linewidth  enter  approximately 
linearly  as  error  terms.  The  curves  indicate  that  very  precise  laser  linewidth 
stability  is  required  to  '-erform  accurate  DIAL  density  measurements  It  is 
expected  that  it  will  be  difficult  to  achieve  pulse-to-pulse  linewidth  and/or 
line  center  frequency  variations  in  the  few  percent  region;  as  a  result,  the 
status  of  laser  technology  should  be  assessed  in  significant  detail  in  this  regard 
before  any  further  commitment  is  made  to  future  development  of  this  density 
measurement  technique. 
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BALLOON  EXPERIMENT  PERFORMANCE 


The  series  of  charCs  in  this  section  provide  data  associated  with  the 
Rayleigh/Mie  Analysis  of  the  balloon-based  experiment.  The  first  two 
charts  on  pages  141  and  143  provide  selected  experiment  parameters  for 
the  500  mJ  and  100  mJ  lasers,  respectively,  which  are  considered  obtainable 
with  existing  off-the-shelf  hardware.  The  baseline  energy,  i.e. ,  the 
500  mJ  and  100  mJ,  is  the  raw  output  energy  of  these  1060  nm  wavelength 
Nd:YAG  lasers  prior  to  wavelength  conversion.  As  a  result  of  converting  the 
1060  nm  wavelength  to  the  353  nm  wavelength,  the  output  energy  is  re¬ 
distributed  to  about  54%  and  9%  of  the  baseline  energy  for  the  1060  nm  and 
353  nm  wavelengths,  respectively.  Backscattering  phase  functions  for  both 
molecules  and  aerosols  at  these  wavelengths  are  provided  in  the  computer 
input  data  shown  in  Appendix  B  on  page  223. 

Receiver  system  parameter  values  such  as  those  shown  for  the  telescope 
diameter,  secondary  mirror  obscuration  area,  filter  characteristics,  and 
PMT  (photomultiplier  tube)  quantum  efficiency  are  the  same  for  both  laser 
configurations;  however,  beam  divergence  and  FOV  parameter  values  were 
varied  to  account  for  eye  safe  criteria  under  day  or  night  operation  con¬ 
ditions.  Although  not  shown,  the  beam  divergence  and  FOV  values  are  the 
same  for  both  configurations  if  the  100  mJ  laser  is  designed  to  also  operate 
at  20  km,  including  the  corresponding  background  interference  values 
(counts/km),  i.e.,  333  counts/km  and  7174  counts/km  for  the  1060  nm  and 
353  nm  wavelengths,  respectively. 
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BALLOON  LIDAR  EXPERIMENT  PARAMETERS 
(500  mJ  BASELINE) 


1060  NM 

353  NM 

fl 

•27J/PULSE 

•045J/PULSE 

PMT  QUANTUM  EFF. 

4% 

30% 

TELESCOPE  TRANSMISSION 

65% 

80% 

TELESCOPE  DIAMETER 

0.5M 

0.5M 

TELESCOPE  AREA  OBSCURATION 

15% 

15% 

FILTER  BANDWIDTH 

0.2NM 

1.0NM 

FILTER  TRANSMISSION 

10% 

7% 

XMTR  BEAM  DIVERGENCE  * 

.6/l.7m  RAD 

.6/l.7m  RAD 

RECEIVER  FOV** 

.8/1. 9m  RAD 

.8/l.9m  RAD 

BACKGROUND (DAY) 

333  COUNTS /KM 

7174  COUNTS/KM 

*  DIVERGENCE  TO  MEET  EYE  SAFE  CRITERIA  FOR  DAY/NIGHT  VIEWING  WITH 
50  MM  BINOCULARS  AND  1060  NM  LASER  OUTPUT  AT  40  KM. 

**  0.2m  RAD  ADDED  TO  BEAM  DIVERGENCE  VALUES  TO  PROVIDE  FOR  MISALIGNMENT 
BETWEEN  TRANSMITTER  AND  RECEIVER  LINE  OF  SIGHT. 
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BALLOON  LIDAR  EXPERIMENT  PARAMETERS 
(100  mJ  BASELINE) 


1060  NM 

353  NM 

el 

.054 J/PULSE 

.009 J/PULSE 

PMT  QUANTUM  EFF. 

4% 

30% 

TELESCOPE  TRANSMISSION 

65% 

80% 

TELESCOPE  DIAMETER 

0.5M 

0.5M 

TELESCOPE  AREA  OBSCURATION 

15% 

15% 

FILTER  BANDWIDTH 

0.2NM 

1.0NM 

FILTER  TRANSMISSION 

10% 

7% 

X  MTR  BEAM  DIVERGENCE* 

.3/. 8m  RAD 

.3/. 8m  RAD 

RECEIVER  FOV** 

.5/1. Om  RAD 

.5/l.0m  RAD 

BACKGROUND  (DAY) 

130  COUNTS /KM 

2802  COUNTS /KM 

*  DIVERGENCE  TO  MEET  EYE  SAFE  CRITERIA  FOR  DAY /NIGHT  VIEWING  WITH 
50  MM  BINOCULARS  AND  1060  NM  LASER  OUTPUT  AT  40  KM. 

**  0.2m  RAD  ADDED  TO  BEAM  DIVERGENCE  VALUES  TO  PROVIDE  FOR  MISALIGNMENT 
BETWEEN  TRANSMITTER  AND  RECEIVER  LINE  OF  SIGHT. 


EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
(FIRING  HORIZONTALLY) 

The  curve  data  on  pages  145,  146,  and  147  provide  expected  errors  In  the 
two  color  density  measurement  technique  at  three  separate  altitudes,  with 
the  lasers  firing  horizontally,  and  no  background  interference.  Single  and 
multiple  shot  noise  effects  and  variations  in  the  uncertainty  in  the  back- 
scatter  ratio  are  also  included  in  the  data.  The  curves  do  show  that 

multiple  shot  averaging  significantly  improves  measurement  accuracy,  in 

— 

particular  in  the  higher  altitude  regions  where  the  (2R  x  /OC  )X 
(1  +  252  contri-bution  t0  error  term  becomes  a  negligible  quantity 

(see  page  101  for  the  error  equation  which  was  used  to  determine  the  error 
in  density  for  the  analysis  presented  in  this  section). 

The  ordinate  values  of  all  solid  curves  in  the  vicinity  of  zero  range  are 
predominantly  due  to  the  errors  associated  with  the  uncertainty  in  the  aerosol 
to  molecular  backscatter  ratio  fiA/jto  353.  The  contribution  of  this  parameter 
to  the  error  equation  remains  essentially  constant  for  the  ranges  shown 
in  these  figures,  and  reaches  its  maximum  value  at  about  the  18  km  altitude 
region,  due  to  an  aerosol  layer  within  the  LOWTRAN  model  (see  page  155  for 
a  better  picture  as  to  the  impact  of  this  parameter  over  the  altitude  regions 
of  interest  in  this  study) .  The  spread  in  the  curves  on  page  146  for  20  km 
reflect  this  impact. 

The  slopes  of  the  curves  for  1000  shots  reflect  the  impact  of  the  uncer¬ 
tainty  in  the  optical  thickness  as  a  function  ?  range  and,  as  stated  in  the 
first  paragraph,  becomes  negligible  in  the  more  rarified  atmosphere  of  the 
higher  altitude  regions. 

Tie  remaining  major  error  contributor  to  the  error  in  density  in  these 
figures  is  due  to  shot  noise  which  is  reflected  in  the  curve  data  by  the 
angular  spread  between  the  single  shot  and  1000  shot  curves.  This  angular 
spread  increases  with  altitude  due  to  a  corresponding  decrease  in  the  number 
of  returned  photons  at  any  fixed  range.  The  slope  of  the  total  error  in 
density  at  any  given  altitude,  therefore,  reflects  the  combined  effects  of 
diminishing  returned  photons  with  range  and  contributions  from  the  previous 
major  error  terms. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  HORIZONTALLY  AT  IOKM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  HORIZONTALLY  AT  20KM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  HORIZONTALLY  AT  3QKM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
(FIRING  UPWARD  30°,  60°  AND  VERTICALLY  DOWNWARD) 

The  curve  data  on  pages  149,  150  and  151  provide  expected  errors  when  the  laser  is 
fired  upward  at  30°.  Page  153  provides  data  for  a  60°  upward  direction.  The 
primary  advantage  to  firing  both  downward  and  upward  is  that  the  balloon  can  operate 
at  a  lower  altitude  and  obtain  very  accurate  data  from  about  5  km  to  high  altitudes 
(an  example  of  this  is  shown  on  page  150.  In  addition,  the  balloon  can  carry 
much  higher  payloads  to  the  lower  altitude,  e.g.,  a  5  km  reduction  in  altitude 
increases  the  payload  capacity  by  a  factor  of  about  3.  On  the  other  hand,  firing 
upward  does  not  duplicate  the  downward  capability  sought  by  an  operational  spacecraft 
instrument,  although  it  does  provide  a  broader  scientific  data  base  for  verification 
of  this  measurement  technique. 

The  curve  data  on  pages  154  and  155  provide  expected  LIDAR  returns  and  expected 
errors,  respectively,  in  the  two  color  density  measurement  techniques  when  the  lasers 
are  fired  downward  from  40  km.  The  data  on  page  154  shows  the  return  counts  also 
for  the  20  and  30  KM  altitudes.  The  data  on  page  155  shows  the  effects 
of  multiple  shot  measurements  on  the  error  in  density.  The  curves  indicate  that 
there  is  basically  no  difference  in  accuracy  between  the  two  lasers  since  atmospheric 
and  backseat  ter  uncertainties  predominate  the  error  equation  under  multiple  shot 
conditions.  The  contribution  from  each  term  in  the  error  equation  can  be  seen  in 

the  curve  data  on  page  165  for  the  100  mJ  laser  conditions  shown  on  page  155. 

It  should  be  pointed  out  that  the  axis  labeled  "Altitude"  on  these  and  subsequent 
charts  is  the  actual  altitude  of  the  scattering  volume  center  as  measured  from 
sea  level. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENS IT?  MEAsnPEMRMrc 
FIRING  UPWARD  30°  AT  10KM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  UPWARD  30°  AT  20KM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  UPWARD  60°  AT  30KM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRING  DOWNWARD  -  906  AT  40KM  (NO  BACKGROUND) 
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EXPECTED  ERROR  IN  TUO  COLOR  DENSITY  MEASUREMENT 
(BACKGROUND  EFFECTS) 

The  curve  data  on  pages  157,  158,  and  159  provide  expected  errors  when 
background  noise  is  introduced  into  the  measurement  of  density.  The 
background  signal  is  based  on  the  FOV's  indicated  in  these  charts. 

Similar  comments  can  be  made  with  these  charts  as  those  made  with  the 
preceding  charts  under  multiple  shot  measurement  conditions,  i.e., 
accuracies  are  well  within  the  10%  design  goal  established  for  this 
study. 

The  contributions  from  each  term  in  the  error  equation  which  were  used 
to  develop  the  data  on  page  159  can  be  seen  in  the  curve  data  on  page  166 
for  the  100  mJ  laser  firing  downward  from  40  km  and  an  0.8  mrad  F0V. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
BACKGROUND  EFFECTS  ON  SINGLE  SHOT  RETURNS 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
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EXPECTED  ERROR  IN  WO  COLOR  DENSITY  MEASUREMENTS 
SAME  BACKGROUND  FOR  TWO  DIFFERENT  ALTITUDES 


*100  mJ  LASER  AT  20  KM  100  mJ  LASER  AT  40  KM 

FIRING  DOWNWARD  AT  -90°  FIRING  DOWNWARD  AT  -90° 

BACKGROUND  FOR  0.8  mrad  FOV  BACKGROUND  FOR  O.B  mrad  FOV 
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AND  LASER  AT  20  KM) 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
(FIRING  UPWARD  FROM  GROUND) 


The  curve  data  on  pages  161  and  163  provide  comparison  data  for  lasers  stationed 
on  the  ground  and  fired  vertically  upward.  Experiment  parameters  used  to  gener¬ 
ate  the  data  shown  on  page  161  are  the  same  as  those  shown  on  pages  141  and  143 
without  any  background  interference;  however,  for  the  data  shown  on  page  163, 
a  higher  energy  laser  (2J  baseline)  and  a  larger  diameter  (1M)  telescope  primary 
were  considered;  all  other  parameters  are  the  same.  This  includes  wavelength 
conversion  efficiencies  (see  page  140),  overall  system  efficiencies  (see  page  223), 
and  backscattering  phase  functions  (see  page  223)  for  the  1060  nm  and  353  nm 
wavelengths.  The  primary  reason  for  the  poor  performance  for  a  ground  based  laser 
is  the  impact  of  the  (2R0t^p  )  X  (1  +  )  term  in  the  error  equation,  and 

which  is  independent  of  the  laser  energy.  The  contribution  of  this  term  can  be 
seen  in  the  curve  data  shown  on  page  167. 


It  should  he  pointed  out  that  these  curves  are  primarily  representative  of  the 
uncertainties  surrounding  the  density  determination  error  equation  and,  therefore, 
do  not  reflect  the  impact  of  additional  information  which  might  be  available  to 
establish  boundary  conditions,  e.g. ,  ground  based  measurements  of  temperature 
and  pressure  can  be  used  to  determine  the  boundary  conditions  for  the  density 
measurements  and,  therefore,  the  expected  error  curves  will  shift  significantly 
downward.  Also,  a  possible  option  for  the  ground  based  instrument  might  be  the 
use  of  a  different  combination  of  wavelengths,  e.g.,  1060  nm  and  530  nm,  which 
should  provide  greater  density  measurement  accuracies,  since  the  impact  of 
uncertainty  (^$^P  )  X  optical  thickness  (2R«)  in  the  error  term  is  less  at 
530  nm  than  at  353  nm.  However,  the  uncertainty  in  the  backscatter  ratio  will 
still  be  a  major  contributor  to  the  error  term.  More  work  is  recommended  in 
this  area  in  order  to  reduce  the  uncertainties  associated  with  this  parameter. 
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EXPECTED  ERROR  IN  TWO  COLOR  DENSITY  MEASUREMENTS 
FIRIM3  UPWARED  90°  FROM  THE  GROUND  (NO  BAQCGROUND ) 
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ERROR  PARTITIONING  IN  ESTIMATED  ERROR  OF  TOO  COLOR 
DENSITY  ME  AS  UK  E.’fENT  TECHNIQUE 

The  following  three  charts  provide  the  individual  errors  expected  from  each 
of  the  contributors  to  the  error  equation  which  was  presented  in  Section  5 
on  page  101.  These  charts  provide  comparison  of  expected  error  data  for  the 
109  n.J  laser  firing  downward  from  40  Um,  with  and  w/o  a  day  background, 
and  firing  upward  from  the  ground. 

Examination  of  the  first  two  charts  indicates  that  the  influence  of  the  term 
containing  the  50%  uncertainty  in  the  backsratter  ratio  has  a  major  influence 
on  the  expected  total  error.  Day  background  has,  of  course,  a  significant  impact 
on  the  expected  error,  although  averaging  with  a  reasonable  number  of  shots  can 
still  limit  the  error  to  less  than  10%. 

Examination  of  the  final  chart  shows  Lhat  the  term  containing  the  uncertainty  in 
the  extinction  parameter  and  total  optical  thickness  traversed  at  353  nm  is  the 
dominant  contributor  to  the  total  expected  error  at  all  altitudes  below  40  km. 
This  term  can  be  significantly  reduced  by  factoring  in  boundary  conditions  at  the 
ground,  i.e.,  ground  pressure  and  temperature  data.  The  anticipated  result  will 
be  that  the  term  containing  shot  noise  at  353  nm  will  dominate  the  expected 
erro  ■  equat  ion  at  higher  altitudes  (also  greater  LIDAR  ranges),  and  the  term 
containing  the  50%  uncertainty  in  the  backscatter  ratio  will  dominate  at  the 
lowe  r  a  1 1  i  tildes. 
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EXPERIMENT  HARDWARE  DESCRIPTIONS 


EXPERIMENT  HARDWARE  DESCRIPTIONS 


This  section  provides  hardware  descriptions  for  the  implementation  of  the  density 
measurement  techniques  under  consideration  in  this  study,  particularly,  for  the 
Rayleigh/Mie  and  O2-DIAL  techniques.  Hardware  descriptions  are  provided  for 
the  Atmospheric  Lidar  Multi-User  Instrument  System  (ALMIS) ,  the  Standard  Test 
Rack  (STR)  LIDAR  Experiment,  as  an  adjunct  to  the  WIND SAT  LIDAR  Experiment,  and 
as  a  proof-of-concept  balloon-based  LIDAR  experiment.  Experiment  hardware  descriptions 
are  discussed  in  the  following  four  subsections  which  provide  pictorial  arrangements 
and  weight,  power,  and  volume  information  for  each  concept.  The  first  three 
of  these  four  hardware  concepts  were  selected  as  viable  alternatives  for  an 
operational  experiment;  although  not  intended  for  operational  use,  the  balloon 
concept  can  provide  a  low  cost  approach  for  a  feasibility  demonstration  experiment 
and,  therefore,  be  a  logical  first  step  towards  obtaining  the  capability  to  develop 
an  experiment  which  can  routinely  provide  accurate  density  measurements  from 
space  on  a  global  scale. 
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WEIGHT  INCREMENTS-  FOR 


STRUCTURE  (INCLUDES  BRIDGE  1480  lb.  (850  lb.  w/o  bridge) 

AND  ALL  PANELS) 

POWER  REGULATOR  /  CONTROLLER  130  lb. 

THERMAL  (HEAT  PIPES,  MLI,  ETC.)  2201b. 

1830  lb.  (1200  lb.  w/o  bridge) 


ACCOMMODATIONS 


50  hours  of  instrument  operation  during  this  time  period.  A  shorter  duty  cycle  will  reduce 
this  requirement  to  values  appropriate  with  a  minimum  number  of  battery  power  modules 
(a  1/12  duty  cycle,  or  5  min.  operation  each  hour,  at  3  KW  requires  only  a  single  power 
module  and  equates  to  an  orbit  average  power  of  only  250  W  which  can  easily  be  handled  by 
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MUST  BE  SHARED  WITH  OTHER  EXPERIMENTS 

CHARGED  TO  PAYLOAD  (SIGNIFICANT  WEIGHT  AND  VOLUME  PENALTY  - 

1632  IbsfKITI 
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•:>!  DESIGN  CHARACTERISTICS  SUMMARY 


Single  STR  C&PH  and  Attitude  Reference  modules  are  also  included  in  the  accompanying  tabulation  and  provide 
complete  autonomy  for  this  experiment.  It  is  conceivable  that  the  Attitude  Reference  module  may  not  be  required 
since  Shuttle  pointing  accuracy  may  be  sufficient,  however,  the  computer  within  this  module,  or  a  compatible 
dedicated  experiment  computer,  should  be  included  with  the  experiment  components. 
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ESTIMATED  POWER  AND  WEIGHT  OF  BASIC  ELECTRONIC 
COMPONENTS  FOR  BASELINE  UINDSAT  CONFIGURATIONS 
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altitude,  i.e.,  the  peak  operating  altitude  would  be  lowered  by  about  0.3  km  with  an  additonal 
100  lb.  battery  pair  and  associated  support  structure  weight. 


ATMOSPHERIC  DENSITIES  FROM  SHUTTLE -BASED  RAMAN  LIDAR  * 

MEASUREMENTS 


SUMMARY: 

Of  the  various  potential  methods  for  measuring  atmospheric  density  from  space 
Shuttle-based  lidar  systems,  this  analysis  considers  only  the  use  of  both 
vibrational  and  rotational  Raman-scattered  return  signals. 

Briefly,  neither  vibrational  nor  rotational  Raman  lidar  return  signals  from  a 
reasonable  lidar  system  appear  to  be  capable  of  good  enough  single  shot  S/N  to 
provide  the  desired  10%  accuracy  at  any  altitude.  However,  if  data  integration 
could  be  satisfactory  for  70  shots,  then  the  calculated  best  night-time  rotational 
Raman  system  (with  height  resolution  of  1  km)  yields  S/N^j20  for  altitude  (Z)^-5  ~n, 
decreasing  to  S/N  =  15  at  Z  =  10km,  S/N  =  7.6  at  Z  =  20  km,  and  S/N  =  3.6  at 
Z  =  30km.  Equivalent  optimum  vibrational  Raman  S/N  values  are  2  to  3  times 
smaller.  An  accuracy  of  10%  would  require  S/N  ^  10. 

Daytime  measurements  are  not  practical  due  to  high  background  levels  yielding 
very  low  S/N  values  (  4^  1).  Eye  safety  considerations  have  not  been 
incorporated  into  these  calculated  results. 


*  General  Electric  PIR  U-1255-IR5-260,  Bethke,  G.W.,  10  November  1978 
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EQUATIONS 
Lidar  Equations 

The  expression  which  relates  the  lidar  receiver  output  to  the  Ildar  system  param¬ 
eters  and  atmospheric  conditions  is  called  the  range  equation.  When  expressed  as 
number  of  backscattered  photons  detected/sec (Ng) ,  the  range  equation  is 


Ns  =  Cn/(8h)]  •  k ( 1 80)  •  As  EL(d/r)2  ET^  (1) 

whcrc  E  =  E0  ;  Epm  (2) 

r  =  ct/2.  (3) 

For  same-wave length  (Mie  plus  Rayleigh)  scattering  (Ag  =  A£) 

k(180)  -  kr  (180)  +  km  (180)  ?  N  •  Qr(180)  +  H  *  QJ180)  (4) 

Ts  *  T)  =  exp  [-/£  (kei  +  kal)  dr]  (5) 

kei  =  ^ri  +  "  NQrj  +  ^ 


Here,  h  is  the  Planck  constant,  k(180)  is  the  atmospheric  backscatter  coefficient. 
As  is  the  wavelength  of  scattered  light,  is  laser  transmitter  output  energy, 
d  is  t'.ic  receiver  objective  diameter,  r  is  the  range  from  receiver  objective  to  the 
propagating  laser  output  energy,  E  is  the  lidar  receiver  efficiency,  T  is  the  one¬ 
way  atmospheric  transmittance,  Eq  is -the  lidar  receiver  optical  efficiency,  Epm 
is  the  detector  quantum  efficiency,  c  is  the  velocity  of  light,  t  is  the  time 
since  lidar  time  zero,  N  is  the  atmospheric  number  density  of  all  gas  molecules, 
Qc(lfO)  is  the  Rayleigh  backscatter  cross  section  for  air,  M  is  the  mass  concen¬ 
tration  of  aerosols,  Q  (180)  is  the  mass- normalized  (Mie)  backscatter  cross  section 
for  aerosols,  kc,  kr,  and  k^,  are  extinction  coefficients,  Qr  is  the  Rayleigh  total 
scatter  cross  section  for  aerosols,  ka  is  the  atmospheric  absorption  coefficient, 
sub-r  and  sub-m  refer  to  Rayleigh  and  Mie  scattering,  respectively,  and  sub-i  and 
sub-r.  refer  to  incident  (laser)  and  scattered  (detected)  wavelengths,  respectively. 
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For  wavelength-shifted  (Raman)  scattering 

k (180)  =  kg(180)  =  Ng  •  Qg(180)  (7) 

=  same  as  above  (eq.  5)  (8) 

Ts  -  exp  [-/J  (kes  +  kas)  dr]  (9) 

kef  =  same  as  above  (eq.  6)  (10) 

^es  a  *Vs  +  ^ms  s  ^rs  +  ^ms  * 


Here,  k„(180)  is  the  Raman  power  backscatter  coefficient  for  gaa  component  g, 

© 

Ng  is  the  number  density  of  atmospheric  gas  component  g,  and  Qg(180)  is  the  Raman 
power  backscatter  cross  section  for  gas  g. 

Daytime  Background  Signal 

Depending  on  the  lidar  pointing  direction,  the  cw  background  signal  will  normally 

be  due  to  diffuse  reflection  of  light  by  the  earth  local  background.  The  number 

of  background  photons  detected  per  second  (n  )  are 

b 

Nb(target)  =  [107  */(16  he)]  •  FsR$Tsxs  B  d2  a2  E  (l2 ) 

for  diffuse  reflection  of  light  incident  on  the  earth  background.  Here  Fs  is  the 
light  flux  (W  ccT^A"*)  incident  on  the  background  at  ^s,  B  is  the  receiver  spectral 
bandpass  (FWHH) ,  c(is  tha  full  angle  field  of  view  of  the  lidar  receiver,  and  Rs 
is  the  background  total  reflectivity  or  local  albedo  at  Xs.  In  the  infrared  at 
3  microns  wavelength,  another  significant  source  of  cw  background  radiation  which 
must  be  considered  is  thermal  emission  from  the  earth's  surface,  air  and  clouds. 

If  Fg  is  the  solar  flux  (actually  irradiance)  outside  the  earth's  atmosphere  and 
Rg  is  the  earth  albedo  from  space,  then  Tg  =1. 
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Lidar  systems  operating  in  the  visble  and  ultraviolet  spectral  regions  usually 
use  photomultiplier  tube  detectors  due  to  their  relatively  high  quantum  efficien¬ 
cies  and  relatively  noise-free  high  gain.  As  a  result,  high  sensitivity  photo¬ 
multiplier  tubes  can  detect  individual  photons  when  illuminated  by  low  light 
levels.  In  this  case,  the  primary  noise  is  due  to  the  statistical  f lunctuations 
in  arrival  of  the  individual  photons (yielding  shot  noise).  Thus  assuming  Poisson 

statistics,  the  standard  deviation  variation  or  noise  count,  N  ,  in  a  specific 

n 


time  (lidar  range)  interval  is 

N  =-/~N  (avg) 


(13) 


If  the  signal  (S)  contains  no  background  photons  and  the  detector  dark  current  is 
negligible  or  discriminated  against,  then 

s/»  *»,/✓»,=  vTT  (14) 

If  the  background  count  level  is  not  negligible,  then 


S/N  -  Ns/Vns  +  Nb  . 


(15) 
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CALCULATIONS 


Using  the  above  equations  plus  the  additional  relations  and  assumptions  discussed 

below  and  sisnmarized  in  Tables  I,  II,  and  III,  Raman  lidar  signals  (N  )  and  signal 

8 

to  noise  ratios  (S/N)  have  been  calculated  for  a  variety  of  potentially  useful 

laser  wavelengths.  Tables  II  and  III  present  optimum  ground  level  (Z*0)  signal 

returns  from  both  vibrational  Raman  and  rotational  Raman  shuttle  -  mounted  lidar 

systems,  respectively.  In  both  cases,  the  single  shot  N  is  tabulated  in  units  of 

s 

photons  detected/micro-second  and  in  photons  detected  (counts)/km  of  range  resolution 
(6.67  micro-seconds).  A  night-time  S/N  is  also  tabulated  for  each  case. 

Table  IV  lists  the  assumptions  and  calculated  results  for  the  best-wavelength 
vibrational  and  rotational  Raman  lidar  systems  (based  on  Table  II  and  III  results) 
as  a  function  of  scatterer  altitude  Z.  Thus  the  vibrational  Raman  lidar  calculation 
is  for  a  near  UV  system  (about  350QA  laser)  while  the  rotational  Raman  lidar  calcu¬ 
lation  is  for  a  mid-visible  system  (5300A  laser).  Here,  the  results  are  tabulated 

as  single  shot  N  /km  of  range  resolution,  plus  night  (and  day)  S/N  for  both  single 
s 

shot  use  and  also  for  70  shots  integrated. 


Raman  Scattering  Cross  Sections 


Although  Raman  cross  sections  of  gases  have  been  measured  at  only  a  very  few 

wavelengths,  the  values  can  be  converted  to  other  wavelengths  as  follows: 

The  Q  (180,  A^)  value  for  desire  incident  wavelength  A-  can  be  calculated  from 
8  / 
a  measured  value  at  wavelength  A-  from  the  relation 

Qg(180,  Xj)  -  Qg(180, 

-Q  (180,  XT)  «  (16) 

8  Ui/Ap-AjfcJ 

where  (*  — >^  )  is  Raman-scattered  frequency  shift  in  cm  *  ,  and  and 

values  are  in  cm. 

In  this  way,  both  the  vibrational  ^  and  the  rotational  ^  Raman  literature  values 
for  Q  (180)  of  nitrogen  and  oxygen  have  been  converted  from  measurements  at  about 

O 

500QA  (see  Table  I)  to  the  wavelengths  of  interest  for  these  calculations  (see 
Tables  II  and  III).  Table  I  also  lists  other  important  information  about  the 
nitrogen  and  oxygen  Ram an  bands.  Finally,  since  the  rotational  Raman  spectra  of 
atmospheric  nitrogen  and  oxygen  (as  well  as  most  other  gases)  occupy  the  same 
wavelength  interval,  the  total  rotational  0^(180)  for  air  is  calculated  as  indicated 
below  Table  1 . 
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Assumptions 

Tables  11,  111,  and  IV  are  calculated  from  equations  1-16  and  the  Information  in 

Table  I  plus  the  following  assumptions  and  literature  values:  The  shuttle  altitude 

is  300  km,  the  laser  transmitter  output  is  1  joule/pulse,  and  the  Ildar  receiver 

2 

is  a  compound  telescope  (2  reflections  plus  1  collimating  lens)  of  1  meter  area 
(d  =  1.128  m.  diam.)  with  an  optical  efficiency  of  0.7.  This  value  of  0.7  does 
not  include  filter  efficiency  Ef  and  any  Raman  band  selection  efficiency.  For  all 
three  tables,  the  filtering  is  assumed  to  consist  of  two  narrow  band  pass  filters 
in  series  (each  2  or  3  cavity  type)  so  as  to  reject  the  relatively  much  stronger 
Rayleigh/Hie  signal.  Thus 

Ef  "  Ef  1  *  Ef  2 

The  detectors  are  all  assumed  to  be  photomultiplier  tubes  with  quantum  efficiency 

E  and  negligible  dark  current  («.l  count/microsecond).  Finally,  the  laser 
pm 

wavelengths  for  were  chosen  on  the  basis  of  existing  or  anticipated  availability 
of  reliable  short-pulsed  high  output  laser  systems. 


Tables  II  and  III  were  calculated  for  only  ground  level  (Z=0)  scattering  and  at  night, 
so  as  to  allow  selection  of  the  best  wavelengths  for  further  consideration.  The 
one-way  Z=0  —  eo  atmospheric  transmittance  values  T^  and  Tg  are  all  based  on  the 

(3) 

very  clear  atmospheric  conditions  of  Allen 

Since  nitrogen  is  the  major  (and  a  well  mixed)  atmospheric  constituent,  the  vibra-. 
tional  Raman  lidar  calculations  of  Table  II  are  based  on  the  nitrogen  Raman  band 
which  has  AV=  2331  cm  ,  Implicit  in  selecting  realistic  Ef  values^  is  the 
assumption  that  the  spectral  filter  bandwidth  B  (FWUM)  is  wide  enough  to  include  the 
entire  nitrogen  vibrat’onal  band.  Thus 

B»300  cm” 1  js  (  V>;)/7 

The  rotational  Raman  lidar  calculations  of  Table  III  are  based  on  the  assumption 
that  the  spectral  filtering  passes  only  the  longest  wavelength  half  of  the  Stokes 
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branch,  with  that  branch  having  601  of  the  total  rotational  Qg(180.  Thus 

Eq  ■  0. 7Ef  *  0. 21Ef  where  EjJ  «  0.6  x  0.5  Ef.  Implicit  in  selecting  Ef  values  is 

the  assumption  that  (FVHM) 

BJS  2  (  y,  (Stokes  max)-  >•) 

As  indicated  earlier,  since  the  rotational  Raman  spectra  of  all  major  atmospheric 
gases  fall  in  the  same  wavelength  interval,  the  rotational  Raman  calculations  are 
based  on  the  total  atmospheric  density  for  N^.and  the  effective  average  rotational 
Qg(180)  (see  beneath  Table  I) . 

The  "best  case"  Raman  lidar  calculations  of  Table  IV  have  the  same  vibrational  Raman 
assumptions  as  Table  II  and  the  same  rotational  Raman  assumptions  as  Table  III, 
except  that  the  atmospheric  density  (N(air))  values  are  here  from  the  U.S.  Standard 
Atmosphere (1962) .  Also  for  Table  IV,  the  range  equation  factor  r  is  a  function  of 
scatterer  altitude  Z. 

r  =  Z(shuttle)  -  Z  (17) 

Since  the  scattering  is  now  calculated  for  various  altitudes,  T.  and  T^  are  functions 
of  Z  which  approach  T^  =  Tg  =  1  at  Z  «=  infinity.  We  assume  that 

k  =  (k  +  k  )  =  C  •  N  (18) 

C  fi 

where  here,  N  =  N(air). 

Thus  T (Z  -  oo  )  =  exp  (-  k  dZ)  =  exp  (-  C  £  N  dZ)  (19) 

where  C  =  -In  T  (0-<»)/(  N  dZ) .  (20) 

Combining  (18)  and  (19)  for  both  T  (Z  -eo)  and  T  (Z  - oo) ,  we  have 

T  i '  Tg  =  exp  [[  £°N  dZ/(  dZ)]  t  [in  Tx  (0- oo)+ln  Ts(0-oo)]  ]■ 

For  Table  IV,  the  T^  •  T^  values  of  equation  (1)  are  calculated  as  a  function  of  Z 
via  equation  (21). 


The  daytime  background  N  and  S/N  numbers  of  Table  IV  are  based  on  a  receiver  field 

b 

of  view  (o<)  of  .001  radian  full  angle.  Since  the  target  irradiance  F'  value  used^^is 

s 

the  solar  irradiance  outside  the  earth's  atmosphere,  and  the  target  reflectivity  R 
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value  used  ^  is  the  earth  albedo,  we  set  T  ■  1  for  the  calculation  of  N.  via 

s  b 

equation  (12). 

CONCLUSIONS 

A  comparison  of  the  results  from  Tables  II  and  III  leads  us  to  conclude  that 
rotational  Raman  lidar  has  the  potential  for  significantly  better  performance  than 
vibrational  Raman  lidar.  Also,  rotational  Raman  lidar  appears  to  perform  best  for 
laser  wavelengths  (  A.  )  in  the  green  and  probably  blue  regions  of  the  visible 
spectrum.  The  best  rotational  Raman  lidar  performance  would  be  found  for  the  near 
UV  (as  with  vibrational  Raman  systems),  except  that  use  of  the  rotational  Raman 
signal  requires  relatively  narrow  very  sharp  cut-off  filters  (or  other  spectral 
isolation)  especially  in  the  UV,  while  band  pass  filter  technology  is  less  advanced 
in  the  UV  as  compared  with  the  visible. 

Table  IV  shows  us  that  for  the  size  lidar  system  and  laser  energy  (1  joule/pulse) 
assumed,  the  "best  case"  green  laser  rotational  Raman  lidar  can  yield  the  necessary 
accuracy  of  10%  (and  thus  S/N >  10  }  during  night  time  only  if  the  1  km  range-resolved 
night  time  signals  are  integrated  over  at  least  16  shots,  and  preferably  more 
(50-100-shots).  During  the  day  time,  the  earth  background  signal  and  noise  dominate, 
making  the  Raman  method  not  useable  during  daytime. 

It  should  be  noted  that  while  the  vibrational  Raman  spectral  requirements  are 
relatively  easy  to  meet  and  thus  the  assumptions  are  reliable,  the  rotational  Raman 
spectral  requirements  are  more  difficult  to  meet  resulting  in  lower  reliability  for 
the  assumed  spectral  purity  and  receiver  efficiency  (Eq) .  Thus  depending  on  the 
exact  technique  used  for  rotational  Raman  spectral  isolation  and  the  "state  of  the 
art"  at  the  time  of  systen  design,  a  rotational  Raman  system  could  perform  signifi¬ 
cantly  better  (it  has  good  potential  in  the  near  UV)  or  worse  than  calculated. 
Finally,  eye  safety  considerations  have  not  been  incorporated  into  these  results, 
and  may  limit  laser  output  to  less  than  1  joule/shot  and/or  require  larger  fields  of 
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view  which  result  in  poorer  filter  performance.  Also,  note  that  the  results  were 
based  on  the  atmospheric  transmittance  for  extremely  clear  air.  With  less  than 
very  clearest  weather,  the  results  for  Z  <  5  km  will  be  degraded,  especially  at 
short  wavelengths. 
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TABLE  I 

RAMAN  BACKS CATTER  CROSS  SECTIONS 
and  related  information 


GAS 

NITROGEN 

OXYGEN 

Vol.  %  in  air 

78.03% 

20.997. 

Stokes  Vibrational  Band: 

Q  (180,  >*)  (cm  2sr  *) 

4.40  x  10"31 

-31 

5.40  x  10  1 

at  (A) 

5145 

5145 

,  -1. 

Band  —  1^.  (cm  ) 

2331 

1557 

Rotational  Bands: 

Stokes  +  Anti  Stokes  total 

.  2  -l 

-29 

,  -29 

Q  (180,  >•)  (cm  sr  ) 
at  >;  (A) 

1.64  x  10 

4.28  x  10 

4880 

4880 

-1 

Rot  Band  Widths*  (FWHM) (cm  ) 

Stokes  Band 

~95 

~  82 

Anti  Stokes  Band 

-  91 

«"*  82 

Band  Max  4V*  (cm  *) 

Stokes  Band 

+58 

+52 

Anti  Stokes  Band 

-56 

-45 

NOTES: 

*  Gas  at  23°C. 

For  Air  Qg(180) 


=  (0.78  x  1.64  +  0.21  x  4.28)  x  10 
=  2.18  x  10  cn>2  sr  ^  at  4880A 
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TABLE  II.  Vibrational  Raman  Single  Shot  Lidar  Returns 
from  ground  level  (Z  =  0)  at  night. 


LASER  SOURCE 

Ruby 

Nd/2 

Nd/3 

Ruby/2 

A;  (A) 

10600 

6943 

5300 

3533 

3472 

K  <A> 

14079 

8284 

6047 

3850 

3778 

Q  (180, N)  (cm2  sr'1) 
& 

1.31 (-32) 

1.09(-31) 

3.84(-31) 

2. 34 (-30) 

2. 52 (-30) 

T. 

l 

0.96 

0.91 

0.81 

0.47 

0.47 

T 

s 

0.98 

0.94 

0.85 

0.58 

0.58 

Ef 

0.15 

0.25 

0.25 

0.1 

0.1 

E 

* 

0.10 

0.15 

0.3 

0.3 

pm 

N^  (counts/ sec) 

— 

0.0221 

0.0686 

0.0843 

0.0891 

Ng  (cts/1  km  At) 

— 

0.148 

0.458 

0.562 

0.594 

S/N  (At  =  1  km) 

— 

0.38 

0.68 

0.75 

0.77 

*  No  photomultiplier  tube  or  Si  photodiode  available  for  this  . 

Assumed: 

Night-time  (no  background) 

E  =  1  ioule/pulse  E  =  0.7  E, 

L  -  r  of 

d  =  1.128  m.  diam  is  for  two  filters  in  series. 

r  =  300  kin 

The  vibrational  Raman  signal  used,  is  that  from  atmospheric  nitrogen,  where: 
AV=V,-^=  2331  cm_1 

N  =  1.951  x  10^  cm  ^  at  20°C.  and  1  atm  (Z=0) . 

8 
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TABLE  III 


Rotational  Raman  Single  Shot  Lidar  Returns 
from  ground  level  (Z  -  0)  at  night 


IASER  SOURCE 

Nd 

Ruby 

Nd/2 

Nd/3 

Ruby/ 2 

h  <»> 

10600 

6943 

5300 

3533 

3472 

Q  (180,  air)  (cm2  s r  *) 

5.79(-3i; 

5. 32 (-30) 

1.57 (-29) 

7. 94(-29) 

8.51(-29) 

T  T 
i  a 

0.96 

0.91 

0.81 

0.47 

0.47 

)^(Stokes  max)-  (A) 

.63. 

27. 

15.8 

7.0 

6.8 

Ef 

.075 

.075 

.06 

.007 

.007 

E 

pm 

.02 

.12 

.20 

.30 

.30 

Ng(counts/yH  sec) 

.0124 

.407 

.961 

.191 

.201 

Ng  (cts/1  km  Ar) 

.083 

2.71 

6.41 

1.27 

1.34 

S/N  (  Ar=l  km) 

.28 

1.65 

2.53 

1.13 

1.16 

Assumed: 

Night  time  (no  background) 

ET  =  1  joule/pulse  E'  “  0.3  E 

L  I  f 

d  =  1.128  m.  diam.  E^  is  for  two  filters  in  series 

r  =  300  km  E  *  0.7EL  =  0.21E 

off 

The  rotational  Raman  signal  is  from  both  atmospheric  nitrogen  and  oxygen, 

where  for  Stokes  plus  anti  Stokes  values, 

Q  (air)  =  0.78  ‘Q  (Nj+0. 21»Q  (0,) 

5  8  *  5  * 

N  =  2.50  x  1019  cm"3at  20°C  and  1  atm  (Z=0) 

8 
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TABLE  IV 


Best  Case  Raman  Lidar  Returns  and  S/N 


VLURAT10NAL  RAMAN 

1 

ROTATIONAL  RAMAN 

z 

(km) 

N(air) 

(cm  3) 

C  NdZ 
C  NdZ 

Ng (cts/km) 

S/N  (night)  1 

N  (cts/km) 

s 

S/N  (night) 

1  shot 

1  shot 

70  shots 

l  shot 

l  shot 

70  shots 

0 

2.55  (19) 

1.. 

.590 

.768 

6.4 

6.56 

2.56 

21.4 

5 

1.53  (19) 

.534 

.671 

.819 

6.8 

4.95 

2.23 

18. *, 

10 

8.60  (18) 

.263 

.555 

.745 

0.2 

3.23 

i  .80 

15.0 

15 

4.05  (IS) 

.121 

.325 

.570 

4.8 

1.67 

1.29 

10. S 

2  > 

1.85  (18) 

.0558 

.  168 

.409 

3.4 

.813 

.902 

7.6 

2  . 

8.33  (17) 

.0251 

.081 

.285 

2.4 

.385 

.620 

5.2 

30 

3.83  (17) 

.0117 

.040 

.199 

1.7 

.  184 

.430 

3.6 

40 

8.31  (16) 

.0032 

.  009 

.097 

mm 

.043 

.208 

1.7 

Day 

earth  N  (cts 

;/0.7/M.s) 

1.88  (6) 

_ 

4.42  (4) 

- 

Da v  S/N  /host  7.)  ! 

— 

.0016 

.013 

— 

.0312 

.26 

ASSUMPTION'S: 

(A) 

3472  to  3533  (avp.  used)  j 

5300 

As 

(A) 

-> 

3778  to  3850  (avp.  used) 

f*  5320 

Q  (180)  (  cm “si 

2.43  (- 

■30) 

1.57  (-: 

29) 

N 

0.78  • 

N(ai r) 

N(air) 

B 

(A) 

50 

16 

E£ 

(2  filter:;  in  r.eries) 

0.  L 

0.2 

K 

o 

0.7  E 

f 

0.21  1 

Ef 

E 

0.  ' 

0.2 

pm 

T  (Z  =••  0+«~) 

0.4/ 

0.81 

Tt  (Z  -•  0 -*'*•) 

0.\S 

0.81 

R 

s 

o.:- 

0.4 

l’s  (  K  on  “A  ^ ) 

, 

i  •  i :  (-■ 

>) 

1.84  (-5) 

Mnn;  Af;.';uu.|> L i i*us : 

(shuttle)  ;  UK >  km  d  -  1.128  in.  iliam. 

r  /(shuttle)  -  7.  At  •  6.6/  s  for  <ir  =  1  I  an 

Ej  I  joti]  i  ■  /  |m  l  d  =  .001  radian  full  angle 

In  ifYinr  dark  current,  (counts)  mu  u>  .Jigilile. 


APPENDIX  B 

SAMPLE  OUTPUT  FROM  COMPUTER  SIMULTATION  PROGRAM 
FOR  RAYLEIGH/MIE  ANALYSIS 


SAMPLE  OUTPUT  FROM  COMPUTER  SIMULATION 
PROGRAM  FOR  RAYLEIGH/MIE  ANALYSIS 
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